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f*"*^ Abstract 

We have applied the mutual inductive method to study higher-order harmonics of complex ac susceptibility x„ = Xn~*Xra 
^~>for YB2C307_5 thin films as function of the temperature and the applied field. The experimental results were compared 
^^ with analytical and numerical results obtained from the Ishida-Mazaki model and the solution of the integral equation 
^^ for the current density, respectively. Both models allow us to reproduced the main experimental features, however, 
.^f- as n increases the numerical model shows notable discrepancies. This failure can be attributed to the current-voltage 
^\) characteristics. Also this investigation yields the activation energy Uc and the critical current density Jc for two samples 

both at T = 0. 
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1. Introduction 

Since the pioneer works of Bean p!] the analysis of the 
harmonics of the complex ac susceptibility Xn = Xn ~ ^Xn 
has become one of the most common methods used for 
the characterization of high-temperature superconductors 
[31 SI . Conventionally, the ac susceptibility technique [5] 
is used in most physics laboratories; however, the induc- 
tive method has shown to be an alternative and highly 
efficient method not only to study the electromagnetic 
properties and the vortex dynamics in superconducting 
thin films but also some important aspects of the models 
[3 El 13 HI IS] • The physical information contained in the 
harmonics of susceptibility have been extensively studied 
[ZllHl Uni HT] . The real part of the fundamental susceptibil- 
ity x' is directly related to the shielding supercurrents (or 
diamagnetic response) flowing in the film; it gives informa- 
tion on how the magnetization in the sample is evolving; 
whereas the energy dissipation occurring in the sample is 
displayed by the imaginary component x" ■ Moreover, it 
is well known that the irreversibility of magnetization is 
the signature for trap flux within the sample. The non- 
linear magnetic response of the superconductor when flux 
pinning is present, is reflected in the presence of higher 
harmonics. Great efforts have been made to understand 
the physics implied in these harmonics not only for super- 
conductors in bulk [ni na nn mi ng im ng but aiso for 
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superconducting thin films [H H^ . The third harmonic X3 
has shown to be an indicative of fiux penetration into the 
sample and has been embraced as a criterion for the de- 
termination of the critical current density Jc El IS] . 
In a similar fashion, the decrease in the value of higher 
harmonics has been related to a reduction of the oxygen 
contents in YB2C307_5 crystals |l|i,. Other works [9^ have 
suggested that models based solely on the critical state 
model [3 120] cannot reproduce the experimental data of 
the higher harmonics due to thermal effects like flux creep. 

One of the motivations of this paper is to test the va- 
lidity of some of the present models on the one hand, and 
determine the values of their parameters on the other. For 
this purpose we aport a study of the higher harmonics 
for the temperature and magnetic held dependence of the 



complex ac susceptibility (xn 



Xn 



«x« 



3, 5) in 



YB2C307_5 thin fllms. The ac field is applied normal to 
the sample surface between and 3500 A/m and the tem- 
perature is varied, for each run, from 70 to 82 K. Then, the 
collected experimental data is compared with the Brandt 
model [m (551 [23] which deals with the integral equation 
for the current density in two dimensions. Furthermore, a 
comparative analysis of the well established Ishida-Mazaki 
model [13j . which is a good approximation for inhomoge- 
neous superconductors, is carried out. This investigation 
also allows us to estimate the activation energy Uc and 
the critical current density Jc both at T = 0. The re- 
sults of the study reveals that the Ishida-Mazaki model 
reproduces better the experimental curves. This suggests 
that the microstructure of thin films behaves as a set of 
layers forming multiconnected network of superconducting 
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Josephson weak links. 

2. Theoretical Models 

In order to interpret our data, we compared them with 
numerical and analytical models. In this article we have 
worked out the numerical calculations derived from the 
solution of the non-linear integral equation for the cur- 
rent density which allows us to analyze the dynamics of 
the critical current density whenever there is a nontriv- 
ial dependence on the position. Also we have adopted 
the analytical expression for the magnetic susceptibility of 
Ishida-Mazaki model which was developed for the study of 
inhomogeneous superconductors. 

In geometrical terms we considered a thin disc of ra- 
dius a and thickness d = 2b immersed in an applied time- 
dependent magnetic field Ba{t) — Sac cos (wi) perpendicu- 
lar to the surface of the film, where Bac is the amplitude of 
the alternating field. These assumptions are valid for both 
models; the main differences are now readily established. 

2.1. Numerical model: current density equation 

Following Brandt [HI im [53], besides the above con- 
siderations, we also assumed a nonlinear response of the 
current-voltage characteristics, i.e., E = p{J)3. The prob- 
lem, then, reduces to find the solution of the integral equa- 
tion for our geometry: 

J{r,t) ^ -^l^' j\\'Q-\ry)[A{r\t) +''-BS)]- (1) 

where A{r,t) = —E[J{r',t)\ and ^o is the magnetic sus- 
ceptibility for the vacuum and the total vector potential is 
given by 

A(r,i) = -Mo r dr' [ dy'Q{v,v')J{y' ,t) -'-Bait) (2) 
Jo Jo ^ 

where r(r, y), r'(r, y). The integrals run over the spatial 
dimensions of the sample and 



Q(r,r') 



-?■' cos (/?' dip' 



[r^ + r'2 — 2rr' cos (/j'Ji/^ 



(3) 



is the integral kernel. 



2.2. Material law 



From the theories of collective creep, flux creep and also 
from the vortex glass picture [Ml [23 [Ml [SB [IHl [IS] , the 
material law equation E{J) becomes strongly nonlinear 
and can be generally described by 



E{J) = Ec exp 



-L/(J) 



knT 



(4) 



where fcs is the Boltzman constant, T is the temperature 
of the sample, the prefactor Ec — BaVo and vq is the flux 



creep velocity. Next we introduce the current density J- 
dependent activation energy 



U{J) 



Uc 



^v- 



(5) 



for the depinning, where Uc is a characteristic energy scale 
and Jc{T, B) is the critical current density which could be 
function of the magnetic field and the temperature. The 
glassy exponent /i is a universal constant in the vortex- 
glass theory with a value ^1. If /i — > then U{J) « 
Uc\^{Jc/ J)- Introducing this into equation Q we obtain 
the power law |22j: 



E{J) - Ec 



(6) 



where m = Uc/{kBT) is a parameter known as the creep 
exponent valid for the interval 1 < m, < oo. Typically 
creep behavior is for the interval 1 < m < 20 [11]. The 
nonlinear resistivity for m ^ 1 is caused by thermally ac- 
tivated depinning with the activation energy ([S]). The case 
m = 1 is the resistive or flux flow regime, which leads to 
normal resistivity (Ohmic behavior) and, to — )■ oo reduces 
to the Bean model. 

The temperature dependence of Jc may be modeled by 
the well known empirical relation [8]: 



Jc(r) = Jco(l-T)^ 



(7) 



where r = T /Tc is the reduced temperature, Tc is the crit- 
ical temperature, Jco is the current density at r = and 
r is a fltting parameter. We also tried other temperature 
dependences for the critical current density, for instance, 
within the context of the Ginzburg-Landau theory the crit- 
ical current is given by Jc{t) = Jco(l — t^)^/^/(H-t^)^/^, 
and for the two fluid model we have Jc{t) — Jco(l— t'^)(1 — 
r^)^". However, we find that the higher harmonic re- 
sponse is well described with equation ([?]). 

On the other hand, the characteristic energy within the 
collective creep model can be described by 



U c{t) ^U coil -T^), 



(8) 



where Ucq is the activation energy at t = 0. With the 
preceding information at hand we could estimate the creep 
parameter. For the temperature range treated in this work 
the creep parameter varies from 1.2 for 82 K to 6 for 70 
K. 

One of the advantages of the power law ^ is that the 
problem reduces to the manipulation of only two parame- 
ters, that is, Uco and Jco- The initial values for Uco and Jco 
were taken from the literature [H [TBI [19] (see next section), 
however, once that the first computation is given, we fixed 
one parameter, say J^o, and manipulate the other. By 
comparison with the experimental results we obtained the 
best values of [/co for each sample. 



2.3. Numerical Integration 

The time integration of equation ^ has to be per- 
formed numericaUy [53]. This may be done by tabulat- 
ing the functions J, E, and A on a 2-D non-equidistant 
grid. A possible nonequidistant grid r^ = {ri,yi) is ob- 
tained by the sustitutions r = r{u) = l/2(3w — u^)a, 
y = y{v) — l/2{3v — v^)b, and then tabulating u = 0, . . . ,1 
and w = 0, . . . , 1 on equidistant grids Uk = {k — 1/2) /Nr 
{k = l,...,Nr) and VI = {l~ Il2)/Ny (I = 1,...,7V^); 
this yields a 2-D grid of A'' = NrNy points with weights 



Wi = WrWy, Wr = drk/dk = 
Wy = dyi/dl = 3/2(1 - vf)blNy, 



3/2(1 - ul)a/Nr, and 
which vanishes at the 
boundaries r — a and y — b. Labeling the points {rk,yi) 
by one index z = 1, 2, . . . , A^ the functions J{r, y, t) become 
time dependent vectors Ji{t) with A^ components and the 
integral kernel becomes a. N x N matrix Qij. Thus, equa- 
tion (fTl) yields the following discrete expression 



N 



Mt) = -f^o^j;T^iQ^jW^) Vi, 



with 



i^j=A,{t)+^Ba{t). 



(9) 



(10) 



This equation can be solved by a Four-order Runge- 
Kutta method or by a step-by-step numerical integra- 
tion; both numerical methods give stable results. In our 
case the space variables were discretized with A^ — 450 
non-equidistant grid points. We used the initial values 
Uco = 480 K and J^o = 3 x lO^A/m^. However, after 
the first computation, we compared the results with the 
experimental ones. We found that the best fitting values 
for UcQ were 930 K and 1300 K for the samples M22 and 
M23, respectively (see next sections). It is to be seen that 
these values were kept fixed during the subsequent com- 
putations, i.e., for the different temperatures. 

Having obtained the sheet current density the total mag- 
netic moment per unit length of a disk was calculated using 
the expression 

/> fa fb 

M^ / d'^rrxJ^27r / r'^dr / dyj{r,y). (11) 
Js Jo Jo 

And, finally, the components of the magnetic susceptibility 
were reckoned using the typical relations: 



-BacTI" JO 



M{t)sm{nujt)diLut), 



BacT^ Jo 



2;r 



M{t) cos{nujt)d{ujt) . 



(12) 



2. A. The Ishida-Mazaki Model 



In the Ishida-Mazaki model [30_ inhomogeneous super- 
conductors are modeled as a stack of layers forming net- 
works of multi-connected loops closed by Josephson weak 



links. In such case the susceptibility is analytically given 

by 



X'l = -(l/47r2)(a-isin2a), 
X'l = (l/47r2)sin2a, 



(13) 



An 



Xn 



1 /sin(n+l)a sin(n — l)a 



Amr^ \ n + 1 n — 1 

1 /cos(n + l)a — 1 



Arnr-^ \ n + 1 

cos(n — l)a — 1 
n- 1 



(14) 



where 



a = 2 sin ^ (sin ( 



^l/2 



sm9 = B^/Ba. (15) 



In this case Bm is the magnetic field which induces super- 
current corresponding to J^ in the loop. This field plays 
the role of the full penetration field H* within the context 
of the Bean model Jj. Also for n =even, x'n ^^'^ Xn do 
not appear. 

On the other hand, we can introduce the temperature 
dependence of the magnetic susceptibility via i?,„. If the 
critical state is present, thus a linear relation between B^ 
and Jc, from equation nh, can be assumed [5]. The param- 
eter r is close to 1 whenever T is more than 10 K below 
Tc and ~ 3.5 for temperatures very close to Tc |9]. For 
the temperature interval of this work the best fitting was 
achieved when r — 1.5, whilst Jco was estimated from the 
experimental curves with r fixed (see section El). 

3. Experimental 

Two superconducting thin films of ¥620307-5 were de- 
posited on 10 mm x 10mm heated LaAlOa substrates, us- 
ing the RF sputtering technique. These films were labeled 
as M22 and M23. Their properties are summarized in Ta- 
ble [H 

Table 1: Properties for the samples used in this work. 



Sample 


Tc 


AT 


Thiekncss 


M70K) 




(K) 


(K) 


(nm) 


xW" A/m^ 


M22 
M23 


86.0 ±0.5 
86.5 ±0.5 


6.0 ±0.5 
5.0 ±0.5 


266 ±5 
260 ± 30 


2.7 ±0.5 
2.5 ±0.5 



For the detection of harmonic generation we have used 
the so-called mutual inductive method [71 [HI Ej ■ In this 
technique a multiturn excitation coil is pressed against the 
film and a pickup coil is coaxially placed on the rear of the 
film to sense the deformed wave. The diameters of the coils 
were approximately 5 mm and the film-coil distance was 
< 1 mm. The driving coil generates a 1 KHz transverse ac 
magnetic field whose amplitude could vary from to 3500 
A/m. A SR830 DSP Stanford lock-in amplifier is used to 



Fourier analyze the deformed signal whose harmonic values 
are fed to a computer program where data are stored and 
displayed in real time. To control the sample temperature, 
films were fixed with silver glue on a Cryomech ST15 cold 
finger. The temperature measurements were performed 
with a 9600-1 Scientific Instruments controller and a Si- 
410NN silicon diode with a resolution of 0.1 K. 



4. Results and discussion 

We have measured the real and imaginary components 
of susceptibility response as function of the applied field 
Ba up to the fifth harmonic for the two superconducting 
films at five different temperatures, namely: 70 K, 75 K, 
78 K, 80 K and 82 K. As mentioned in the introduction of 
this contribution we shall focus on the analysis of harmon- 
ics for n = 3,5. However, in order to verify the consistency 
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Figure 1: Representative curves for the fundamental harmonic of ac 
susceptibility response as function of the Ba at 78 K. (a) experimen- 
tal curves for the two samples used in this work, (b) Ishida-Mazaki 
calculations and Brandt model. There is an overall agreement be- 
tween theory and experimental data. 



of the experimental data among the films and our models 
we have plotted in Fig. IT] representative curves for the 
fundamental component of susceptibility at 78 K. Figure 
[Ik compares the curves for the real and imaginary com- 
ponents for the two samples. We can see that there is a 
minor difference in the structure of the curves. Figure [lb 
compares the curves for the real and imaginary compo- 
nents for the two models used in this work. As the field 



increases the amplitude of the numerical curves for x" falls 
off slower to zero than that of the Ishida-Mazaki calcula- 
tions. Similarly, the opposite occurs for x'- These slight 
discrepancies arise from the different assumptions for each 
model. In particular the numerical calculations assume 
creep effects. The complete analysis for the fundamen- 
tal harmonic has been extensively study by Acosta et al. 

The field dependence for the third harmonic of suscep- 
tibility is shown Fig. [2| It is to be seen that the peaks of 
sample M22 are shifted towards the left for a given temper- 
ature in comparison with M23 [see Fig. [2]a-d] . The physics 
involved in these curves is as follows. If the only mecha- 
nism governing the flux dynamics is pinning then for suffi- 
ciently low temperatures the vortices are strongly pinned 
and the maximum current density is achieved. As tem- 
perature increases vortices star to depin and the current 
density is significantly reduced. For temperatures close to 
Tc the sample goes to the critical state unless the super- 
conducting state is lost. Thus the slight difference in the 
position of the peaks for the two samples is due mainly to 
the lower Tc of M22. This means that for a given field and 
temperature sample M22 goes faster into the mixed stated 
than sample M23. 

Notice also that as the temperature decreases the tran- 
sition in x's and x'l, broads. Moreover, the peak height in 
both components slightly decreases for temperatures far 
from Tc but there is a great augmentation close to T^. 
The increase in the magnitudes of the peaks reflects the 
increase in the nonlinearity in the magnetization, includ- 
ing hysteresis, as a function of the applied field. Another 
contribution to the increase of the peak comes from the 
field dependence of the magnetic susceptibility [see equa- 
tion (12)]. As the amplitude of the applied field tends to 



zero the components of the magnetic susceptibility diverge. 
For comparison with the experimental data we have 
plotted the theoretical results for the third harmonic re- 
sponse for M23 in Figure [2] e-h. The curves for M22 (not 
shown) only differ by a scaling factor that causes a reduc- 
tion in the broadening but the curve structure is of same 
nature as that of M23. It is evident that the analytical 
calculations are in better agreement with the experimen- 
tal results. Based on these results the best fitting values of 
Jco for the Ishida-Mazaki model were 2 x 10^^ A/m^ and 
8 X 10" A/m^ for M22 and M23, respectively. Further- 
more, as can be seen from this figure the amplitude of the 
curves remains constant; this is attributed to the fact that 



the magnetic susceptibility ( 14 ) does not contemplate the 



variation of the amplitude of the applied magnetic field 

On the contrary, the numerical calculations reproduce 
better the amplitude of the experimental curves but the 
peak does not shift proportionally to the temperature 
changes [Fig. [2] g-h] . Also the curves do not fall off as 
fast as the experimental ones. These effects are probably 
due to the simplification realized in the power law relation 
([m which cannot account for the temperature changes in 
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Figure 2: Third harmonic of ac susceptibility response as function of applied field for the whole set of temperatures, (a) and (b): experimental 
curves for sample M22. (c) and (d): experimental curves for sample M23. (e) and (f): Ishida-Mazaki calculations for M23. (g) and (h); 
Brandt numerical calculations for M23. Left panels for real components. Right panels for imaginary components. As temperature lowers the 
main peak shifts to higher fields. The Brandt model qualitatively reproduces the shape of the experimental curves but the peak does not 
shift proportionally to the temperature changes. 
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Figure 3; Fifth harmonic of ac susceptibihty response as function of applied field for the whole set of temperatures, (a) and (b): experimental 
curves for sample M22. (c) and (d); experimental curves for sample M23. (e) and (f) Ishida-Mazaki calculations for M23. And (g) and (h) 
Brandt calculations for M23. Left panels for real components. Right panels for imaginary components. As temperature lowers the main peak 
shifts to higher fields. The Brandt model shows a higher deviation in the agreement with the experimental curves. The peak also does not 
shift proportionally to the temperature changes. 



the same proportion. This may be so, since the model we 
used only evaluates the resistivity as function of the cur- 
rent density which in turn is function of the vortex mobility 
in the film [5T1 [22 [23] . This may indicate that the model 
is not evaluating some parameter which is important to 
reproduce the temperature shift. 

In Figure [3] we plot the fifth component of susceptibility 
for the real and imaginary components for the two sam- 
ples as well as the analytical and numerical calculations as 
function of the temperature and the applied field Ba ■ For 
the same reasons as in the third harmonic we have only 
presented the theoretical curves for M23. 

Here again we can observe that for temperatures close to 
Tc the peak height greatly increases for the sample M22. 
Also as the temperature goes down from Tc the curves 
become broad and the peak height is more stable. For 
comparison, the broadening for the curves of M23 is appre- 
ciably wider than that of M22 and the peak height, as the 
temperature goes down, increases. The plots also exhibit 
an oscillatory behavior for low magnetic fields. The anal- 
ysis of the model has shown that this is mainly caused by 
the abrupt variations of the magnetic susceptibility as the 
applied field Ba goes to zero through the variable a. From 



equation (15) we can clearly evidence that for Ba ^ B,n 
the function a varies very rapidly. Thus we can confirm 
that the overall structure of the data is in excellent agree- 
ment with the analytical model. 

The Brandt model, on the other hand, can only repro- 
duce the overall structure of the experimental curves [see 
Fig. Islg-h]. Whilst, at low magnetic fields, the oscillatory 
behavior is not present. Notice that the curves do not 
shift proportionally to the temperature changes. Also the 
peak height of the real component increases as the tem- 
perature increases. In order to verify the cause of these 
discrepancies we tried several temperature dependences of 
Jc- Similarly, we manipulate the value of the creep pa- 
rameter up to rra = 100, this manipulation is equivalent 
to a variation of the activation energy Uc- By following 
this procedure we are, in fact, inducing the model into 
the flux creep regime and the critical state, nevertheless, 
no new features appeared. These tests also support our 
previous arguments given for the third harmonic that the 
unexplainable effects may imply that the power law (pi) 
for the current-voltage characteristics in the flux creep is 
not valid for the range of temperatures treated in this ex- 
periment. We have seen that the whole structure of even 
harmonics is well reproduced but perhaps a change in the 
material law |19| . or a magnetic field dependence of Jc or 
Uc may take place [121 [SI]. However, in this article we did 
not include this models since there are parameters that 
we cannot measure, like, for instance, the characteristic 
length ^ of the flux lattice line. 

5. Conclusions 

Higher-order harmonics of the complex ac magnetic sus- 
ceptibility have been measured for two superconducting 



thin films. We have compared the experimental results 
with theoretical ones which were based on the Ishida- 
Mazaki and the Brandt models. We have shown that 
the overall behavior of the analytical model is in excel- 
lent agreement with the experimental data. Based on this 
vision we may conclude that our thin films behave more as 
network of superconducting layers coupled by Josephson 
weak links. 

On the other hand, although the overall structure of the 
magnetic behavior was reproduced by the Brandt model, 
it could not explain the peak shift as function of temper- 
ature. We attribute this failure to the simplification that 
we have made of the material law (|4|. Alternative ap- 
proaches such as magnetic field dependence of the critical 
current density and activation energy must be considered 
to account for the lack of agreement with the experimental 
results. 
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